A0-A1M  4(1 
UNCLASSIFIED 


DEMONSTRATION  OF  A  JACORI  POLVNOHIAL  REPRESENTATION  OF  1/1 
A  KRONECKER  DELTA. .  <U>  JOHNS  HOPKINS  UNIV  LAUREL  HD 
APPLIED  PMVSICS  LAD  R  A  FARRELL  ET  RL.  SEP  8S 
JHU/APL-T8-13M  NRM24-R9-C-52D1  F/B  12/1  NL 


MiCrfOCCPV  »E  SOLUTION  tes 


A)D-A/6o  V6o 


.iHU  API 
TG  1350 

SEPTEMBER  1985 


DEMONSTRATION  OF  A 

JACOBI  POLYNOMIAL  REPRESENTATION 

OF  A  KRONECKER  DELTA  FUNCTION 


RICHARD  A.  FARRELL 
ERNEST  P.  GRAY 
ROBERT  W.  HART 


f?v 


A’tj;  }  frit  puhl'C  c.  j  *»  =  ,  1 


B5  10  22  042 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


ia  REPORT  SECURITY  CLASSIFICATION 

l  nda^ified 

2a  SECURITY  CLASSIFICATION  AUTHORITY 

2b  DEC  LASS!  F  iC  AT;  ON  DOWNGRADING  SCHEDULE 

n  a 

4  PERFORMING  ORGANIZATION  NUMBER'S- 

J HI  API  I  < i  13<0 

6a  NAME  OF  PERFORMING  ORGANtZAGON 


PLEASE  FOLD  BACK  IF  NOT  NEEDED 
FOR  BIBLIOGRAPHIC  PURPOSES 


REPORT  DOCUMENTATION  PAGE 

|  lb  RESTRICTIVE  MARKINGS 


j  3  distribution  availability  of  report 

Approved  for  public  release;  distribution  is  unlimited 


S  MONITORING  ORGANIZATION  REPORT  NUMBER'S, 

I  HI  API  TG  1350 


6b  OFFICE  SYMBOL  j  7a  NAME  OF  MONITORING  ORGANIZATION 

i  it  Applicable  > 

1  ^  Nasal  Plant  Representative  Office 

j  7b  ADDRESS  iC‘t y.  S/afe.  ana  ZIP  Code  i 

j  Johns  Hopkins  Road 
I  l  aurel,  Maryland  20707 


f  he  lofms  Hopkins  l  rmer-itv 
Applied  I’hvsus  l  ahot.itorv 

6<  ADDRESS  ittate  and  C:P  Cole  ■ 

lohns  Hopkim  Road 
l  aurel.  Mars  land  2<n>~ 

na  NAME  of  c oNDING  SPONSORING  ORGANIZATION 

Naval  Sea  Ss stems  Command 

6'.  ADDRESS  C  *►  $?a.'e  and  Zip  Code . 


U  ashmt'ton.  IX  20362  | 

M  ud''  Secuftt*  Classification , 

Demonstration  ol  a  Jacobi  Polynomial  Representation  of  a  kronecker  Delta  Junction  (U) 

12  PERSONAL  AUrHQRl5' 

larrell.  Richard  A  .  Gray,  Ernest  P  ;  Hart,  Robert  W. 
i  i  :~a  t*p?  0C  RfP('RT  "  "  Tl3b  TiMECOvERED  [  14  DATE  OF  Rf 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

N00024-85-C-5301 

■ 

L 

10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

' 

_  J 

NO 

NO 

X83J 

_ 

ACCESSION  NO 

I  echnnal  Memorandum 

, ,J  Pf * f  M  l  N  ’  A  n  V  N o  r  A  rrQN 


1  K( >S1 


14  DATE  OF  REPORT  (Year  Month  Day )  16  PAGE  COUNT 

1985,  September  14 


:  Os A T i  CODES 


,  18  SUBJECT  TERMS 


vjRt-,,-p  +  subgroup  _ j  Kronecker  delta  function 


,  j  Jacobi  polynomials 

— -  -  - — 'j  fluid  shells  ^ 

•  ••  •*•*»>■•.*  '  necessary  and  identify  by  bfock  number  i 

\  1  i  derr.ed  In  Noble  has  been  generalized  to  obtain  an  expansion  of  the  Kronecker  delta  function  as  an  infinite  series  involving  the  products 

■  •  l(i  noimals 


’RA-;  :  I  21  ABSTRACT  SECURITY  CLASSIFICATION 

□  same  as  rpt  □  dtic  USERS  I  Unclassified 

j  22b  TELEPHONE  < Include  Area  Code ) 

(301)  953-5000 


22c  OFFICE  SYMBOL 

NPRO 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


JHU/APL 
TG  1350 

SEPTEMBER  1985 


Technical  Memorandum 

DEMONSTRATION  OF  A 

JACOBI  POLYNOMIAL  REPRESENTATION 

OF  A  KRONECKER  DELTA  FUNCTION 

RICHARD  A.  FARRELL 
ERNEST  P.  GRAY 
ROBERT  W.  HART 


:  THE  JOHNS  HOPKINS  UNIVERSITY  ■  APPLIED  PHYSICS  LABORATORY 

j  Johns  Hopkins  Road.  Laurel.  Maryland  20707 

. . 1  8‘  i'  '>3Cl  A  I*'  T »»♦.-  D»‘; •.»*»*  -  *  *hp  Na.y 

I 

I 

I  Approved  for  public  release,  distribution  is  unlimited 


:i>n\S  fMPM\S  UNI VI  HSITV 


APPLIED  PHYSICS  LABORATORY 

;  A,  Hti  VARN  I  ANP 


ABSTRACT 


\  formula  derived  by  Noble  has  been  generalized  to  obtain  an  expansion 
ol  the  Kroneeker  delta  function  as  an  infinite  series  involving  the  products  ot 
two  Jacobi  polynomials. 
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l AUBLL  MABVL AND 


INTRODUCTION 


In  our  investigations  of  the  inviscid,  incompressible  How  of  rotating  fluid 
shells  confined  between  concentric,  spherical,  rigid,  co-rotating  boundaries, 
we  encountered  the  need  for  a  proof  of  a  Kronecker  delta  function  representa¬ 
tion  in  terms  of  a  series  of  products  of  two  Jacobi  polynomials.  In  particular, 
we  required  a  proof  that 


bi.i  i 


(/*+ 

'  2“  '  " 


(A)  , 


(1) 


w  here 


.  i.u  -  /( A) 


£ 


*  a  (I 


(*  +  Ji/2  +  //2+  '/i ) 

(k  +  p  +  /+  Vi)(k  +  n+ 


x  P^'1,  •  "(\)  P'Sj1'  ■  "  (A)  ,  0  <  IAI  <1. 


(2) 


For  («,d)>  1,  the  A)  are  the  Jacobi  polynomials  which  are  orthogo¬ 

nal  over  the  interval  (-1.+  1)  with  the  integration  weight  factor  w(x)  = 
(1  —  a )" ( 1  +.V)’1.  They  are  normalized  by  the  relation  P„u'^(\)  =  T (a  +  n+  1)/ 
[Hn  +  1  )T(a  +  1)],  with  T(z)  the  gamma  (factorial)  function.  For  ( a,0 )  <  -  1 , 
the  Jacobi  polynomials1  P„(<,",)( A)  are  defined  in  terms  of  the  polynomials  for 
(n.li)  >  -  1  through  repeated  applications  of  the  contiguous  relations1 ' 

(2 n  +  n+ii)  P„‘"  u'  (A)  =  ( n  +  n  +  ii )  P.,1"'”  (A)  -  (n  +  tf)  P'"'f  (A)  (3a) 


and 


(2n  +  „  +  ,i)  P,;  "J  "(A)  =  (n  +  «r(i)P'"J'  +  (/» +  «)Pjr|’ ( A)  (3b) 
with  P,!"'"’  (A'  -  1  and  P'Y'(A)  =  0. 

Certain  sums  involving  products  of  two  Jacobi  polynomials  have  been  previ¬ 
ously  evaluated  (see,  for  example,  Refs.  3  -  6).  Our  method  for  evaluating  the 
sum  in  Eq.  2  employs  mathematical  procedures  analogous  to  those  which 
Noble'  used  to  show  that6 
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2‘  '  '(1  -  v)  ‘'(l+.v)  V-v)"^  ‘  'A(x.y) 


-  *!r(*  +  c+l)(2*  +  c+l) 

r (a + b~ r)  - 7— — 

r(Ar4"<?+l)r(Ar  +  Zj+l) 


pku  h-h>(x)  P*,a‘  "’O’)  ,  (4) 


where 


I  1  if  -  1  <  y  <  x  <  1 

A(.y,v)  =  .  (5) 

I  0  if  -  1  <  x  <  v  <  1  . 


PROOF 


The  proof  is  done  in  two  stages:  first  the  /  =  0  case  is  treated,  and  then 
the  /  *  0  case  is  analyzed.  Both  stages  use  the  expansion  of  a  function  of  two 
variables  in  an  infinite  series  of  Jacobi  polynomials  of  one  variable  and  coeffi¬ 
cients  that  depend  on  the  other  variable.  In  particular,  the  coefficients  in  the 
expansion  are  obtained  by  respresenting  the  Jacobi  polynomial  as  a  finite  se¬ 
ries,  evaluating  the  resulting  elementary  integrals,  and  re-summing  the  resul¬ 
tant  series. 

In  proving  Eq.  1,  it  is  convenient  to  use  the  relationship 

p 'n'"  V)  =  (l+*>  for  n  >  1,  (6) 

in  order  to  write  Eq.  2  as 


1 


H  ♦  I, ft  -  I 


(2n+  1) 


P, 


<*.-/+ 1 it 


(A) 


+  LLTA!  y'  Vi )  ()l+/4. 1  ..I, 

k(k  +  l) 


8 


*  ,  (A)  PkYi-\  :'l) (A)  , 


(7) 


where  the  first  term  is  the  A:  =  0  term  of  the  sum  appearing  in  Eq.  2. 
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CASE  OF  /  =  0 

For  /  =  0,  we  expand 


Hw(x,y) 


- -  l  (1  -/) -("  +  3,2V/ 

(1  +  a)  ( 1  +y)  J  -I 

£  Cti0'(y)  /V“+  -"(a), 

*=  0 


where 


(8a) 


vv 


x  if  y  >  x 
y  if  y  s  x  . 


(8b) 


Multiplying  Eq.  8a  by  (1  +  a)(1  -a)"*  ''  P„**""u(x)  and  integrating  from  a 
=  -  I  to  1  yields 


(m+l)2"  +  ?,-Cn(>’) 


)||*  V2 


( m  +  n  +  3/2 )  (2m  +  n  +  5/2 )  (I+.V) 


dx(l  ~x)“*  :  '"(a)  dt(\-t) 


\: 


■  (n  +  3/2) 


,*.1 


dx{  I  -  Vp*  P, 


lu*  ;.l) 


(AT) 


(1-/) 


=  (  I  +  V)  '{fl+'A) 


dx(\-x)“*'!  P, 


lp+  ■'.!) 


(A) 


+  2“ 


■“i: 


dxP**':M(x) 


+  :(l-y) 


ip 


.  \  P  ♦  ’  -  p  < p  +  1 : ,  I  ) 


dx(\-x)“+  P, 


(A) 


(9) 


The  integrals  are  readily  evaluated  by  expanding  the  Jacobi  polynomial  as1 


m 

P m (a)  =  Y)  y(m,a,(3;r)(\ -x)r. 


(10a) 
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with 


y(m,ct,f3;r)  =  r(a  +  m  +  1  )T(a  +  0  +  m  +  r+  1 )  (  -  Vi )'/ 

[/•!(m-r)!r(a  +  r  + 1)  T(a  +&  +  m+  1 )];  (10b) 

integrating  each  term;  and  re-summing  the  series.  In  particular. 


(m- t- 1^**5'2 

( m  +  ^  +  3/2 )  ( 2m  +  /i  +  5 /2 ) 


C\Td)  =  (1  +y)-'(li+'/i)-' 


x  E  +  1/2 •  * .'' ) 


r 


-2 


2^  +  5/2 


[d-^)r+,-2r+l] 


c 


(r  +  n  +  3/2)  (r+  1) 


+  2 


,<2 


2r  +  “  +  7/2-2‘,  + 5/2(  1  -_y) 


r+  I 


(H+'/2) 


(r  +  1 )  (r  +  /x  +  3/2) 


=  —  2M 7/2 ( 1  +  j)  ~l  (m  +  n  +  3/2)  “2 


- - 1  ’ ( - 1 )  -  PL\T’-l'(y) 


2<i  +  5/2 


(w+  \  )(m  +  n  +  V2) 


P  (p+  '/!.!) 

*  m 


O')  . 


where  the  fact  that 


01) 


7(m  +  1  ,fi  +  Vi,  -  1 ;r  +  1 )  = 

-  (m  +  n  +  3/2)27  ( +  Vi , 1  ;r)/[2(r  +  p  +  3/2)(r  +  1 )] 

is  used  to  arrive  at  the  next  to  the  last  line,  and  Eq.  6  is  used  in  the  final  step. 
Inserting  this  result  in  Eq.  8a  yields 
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H{0)  (x,y) 


2**  +  5/2 

(l+y)(l+x) 


(1  -tv) 

( /i  +  Vi ) 


-  e  p^"-u(x) 

1=0  (^+  1  ) 


l  /  2  y*_;  I  t  ( I  +y) ( I  +-*■) 
(2/i  +  1 )  \  1  —  vv/  “  (2/1+1)  8 


OD 

x  E 


(2/t  +  /l+  ’/2) 


p  (  m  ♦  1  • .  u 

‘k  -  I 


<»  W'1  ’  U) , 


(i: 


(i; 


with  tv  defined  by  Eq.  8b.  Since  Poa'ff)  (A)  =  1,  comparing  Eq.  7  with  / 
0  to  Eq.  12b  with  jr=^  =  A  yields 


1 

(2/t  +  l) 


( 


which  demonstrates  that  Eq.  1  holds  for  /  =  0. 


CASE  /  *  0 

For  I  ^  0,  we  consider 


//(/>  (jc,^)  = 


'  0  for  1  >  y  >  x  >  -  1, 

r(/i  +  5/2) 


r(/i-/+5/2)/!(i  +x)  ay 


d'-i 

— |  [(1  -J')  ~<"+3/2>(Jc-^)/] 


for  1  >*>>’>  - 1  . 


Application  of  Leibnitz’s  rule  for  differentiating  a  product  yields 


.i+i 


//(,W)  =  is, (x.  y\  - 1/<^) . 


/  1  4. 
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with  S,  = 


(M+3/2)r(M+/+3/2)r(/)(-n,-'-,(i->-)-<''+3/2+'> 


r(/i-/+5/2>r{/+2)r(/+  i)r(/-/) 


and  U(x,y)  =  1  if  1  >  x  >  y  >  —  1,  and  U(x,y)  —  0  otherwise.  The  func¬ 
tion  HU)(x,y )  is  expanded  in  terms  of  Jacobi  polynomials 


HU)(x,y )  =  £  CK0)(y) 


(x) 


(16) 


k =0 


with  the  coefficients  determined  from 


(m  +  1)2m“/+5/2 


Cin(y) 


(m  +  n  —  l+3/2)  (2m  +  1+5/2) 

(1+j <)P^-,+  v'-"(x)  dx 

E  s>  M-/+W.l;r)  j  . 


|  HU)(x,y)  ( 1  —  1/1 '  *  -* -v\p(«-/+ '4.)), 


/-I  m 

E  s<  I 

(=0  r=0 


(17) 


The  integral  is  evaluated  by  changing  the  variable  from  x  to  z  *  (x-^)  /  ( 1  -y), 
which  leads  to 


I(m.l)  =  E  y(m,n~l+  '/2,l;r)(l  ->>)"~'+w+,+2r(M-/+r+3/2) 
T(f  +  2)S(  (!->-)' 


x  E 


1  =  0 


Y(n—l  +  r  +  t  +  l/2) 


(18) 


Substituting  for  S, ,  the  inner  summation  appearing  in  Eq.  18  reduces  to 


(#1  +  3/2)(-l)'-1  l  T(M  +  3/2  +  /) 

T(/i  — /+5/2)  ,r0  /  rU  +  3/2  +  /-(/-2-r)J 


x 


(/-!)! 

/!  (/—  1  -t)\ 


(-1)'  • 


(19) 


For  /  >  r  +  2,  we  note  that 
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o(r,l)  = 


(M  +  3/2H-1)-1  -z)1-')  =  0  for  r  >  0. 

r(n-/+5/2)  bz'~2~r  '*-■  (i 


The  derivative  is  zero  at  z  =  1  because  (/-  2-r)  <  /  -  1.  The  fact  that  the 
derivative  reproduces  the  series  in  Eq.  19  can  be  readily  verified  by  writing 
the  factor  (1  -z)‘  1  as  a  polynomial  in  z,  differentiating  each  term,  and  then 
setting  z=  1.  For  /  <  r  +  2,  we  note  that 


a(r,l)  = 


(~1)/~i(m  +  3/2) 
r(M-/+5/2) 


f*  p  p-i 

\  ds,  \  ds2.  .  .  \  dsjSj 
Jo  Jo  Jo 


/,+'1  (1 


(-i)'-lr(n+s/2)r(r+i) 
r(n-l+5/2)r(r  +  2-l)T(n  +  5/2  +  r)  ' 


The  multiple  integral  in  Eq.  21  is  readily  evaluated  by  recalling  that  integra¬ 
tion  by  parts  leads  to 


I.*'  L'*!  -  *  1  ,<fc' - 


The  fact  that  the  multiple  integral  reproduces  the  series  in  Eq.  19  can  be 
readily  verified  by  writing  (1  -s,  ),_1  as  a  polynomial  in  s,  and  integrating 
each  term  separately.  Inserting  Eqs.  20  and  21  into  Eq.  18  yields 


(-l)/-'r(M  +  5/2)  " 

IthTmT 


r  (n~l+r+  3/2)r! 
x  ro+s/z+uir-M -01 


otherwise. 


For  m  >  /  -  1,  using  the  definition  of  y(m,n~l+  Vi ,  1  ;r)  and  introducing 
the  indices  r’  =r-(/-  1)  and  m’  =m-(l-  1)  leads  to 
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I(m,l)  = 


(m'  +n+  Vi)  rOt  +  5/2-/) 


X 


m  ‘ 


£ 


[(j'-l)/2]rT(m'+M  +  /+l/2+r<) 
(m'-r')!r'!r(/i  +  /+3/2  +  r') 


_ r(^  +  5/2)  ( v^)7-1 _ 

(m  +  ft+3/2)(m  +  /i-/+3/2)r(n-/+5/2) 


P  ( \i  +  /  +  ’/l.  —  I  ) 

*  m  -  ( /  -  1 ) 


O')  . 


(24) 


where  the  definition  of  y(m-l+  1,  n  +  l+  Vi,  -  1  ;r)  was  employed  to  obtain 
the  last  line.  Inserting  Eqs.  23  and  24  into  Eqs.  17  and  16  yields 


2-ili+m)r(n  +  5/2)  “  (2m  +  fi-/+5/2) 

T (n  +  5/2  —  l)  mtf_t  (m+l)(m  +  M  +  3/2) 

x  pr('^Vrnu)  p» .  (25) 


Changing  the  summation  index  to  m '  =  m  -  (/- 1 ),  writing  the  w'  =0  term 
separately,  and  using  Eq.  6  in  both  the  m'  =0  and  the  m'  ^ 0  terms  lead  to 


2-(,‘-'/,,r(*t+5/2)  r o+a-)(i+^) 
(l+x)r(M  +  5/2-/)  L  (2^+T)  +  8 


X 


(2  m’  +l+n+  Vi) 
m’(m'  +  /) 


p  ( n  +  /  +  Vi ,  I )  /  V  p  ( p  -  /  +  ‘/J ,  1 ) 

"  m’-i  '  +  /  -  1 


■  (26) 


Setting  x=y  =  \  in  Eq.  26  and  comparing  with  Eq.  7,  one  finds 


H{,)(  A,  A) 


2-(M-'/.,r(M  +  5/2) 
d  +  A)r(^  +  5/2-/) 


0,  /*()  , 


(27) 


where  the  vanishing  of  HU)  ( A,  A)  follows  directly  from  the  form  of  Eq.  15. 
This  completes  the  proof  of  Eq.  1. 
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